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Abstract 

We calculate the local density of states (LDOS) of a superconductor-normal metal sandwich at arbitrary impurity 
concentration. The presence of the superconductor induces a gap in the normal metal spectrum that is proportional 
to the inverse of the elastic mean free path I for rather clean systems. For a mean free path much shorter than the 
thickness of the normal metal, we find a gap size proportional to / that approaches the behavior predicted by the 
Usadel equation (diffusive limit). 
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A superconductor that is in good metallic con- 
tact to a normal metal induces a finite pair ampli- 
tude on the normal side. As a result, the normal 
metal acquires superconducting properties like in- 
finite conductance and the Meissner effect (see [1] 
and references therein) . The presence of the super- 
conductor also has consequences on the spectrum: 
the constant density of states around the Fermi en- 
ergy of a bulk normal metal is strongly modified 
and shows a pseudogap for a clean normal metal. 
That is, it vanishes at the Fermi energy and rises 
linearly close to it [2]. Another well-known result 
on the spectrum has been obtained in the dirty 
(diffusive) limit [3,4]: in this case, there is a gap 
in the spectrum (called minigap to avoid confusion 
with the gap of the superconductor) that has the 
same order of magnitude as the Thouless energy 
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Erh = hD/d 2 , here D is the diffusion constant of 
the normal metal and d its thickness. 

In related work, the spectrum of a ballistic 
normal cavity connected to a superconductor has 
been studied. A classically integrable cavity shows 
a LDOS that is linear in energy, whereas a chaotic 
cavity exhibits a gap in the spectrum [5] . 

The goal of the present work is to investigate 
a (moderately) disordered normal metal and to 
study in detail how the linear rise of the LDOS for 
the clean system transforms into the minigap in 
the diffusive system. To this end, we have solved 
the real-time Eilenberger equation for the quasi- 
classical 2x2 matrix Green's function g numer- 
ically for the planar geometry shown in the in- 
set of Fig. 2 (see [1] for additional details of this 
method). Impurity scattering is taken into account 
by including an impurity self-energy of the form 
(g) /2r (Born approximation) determined in a self- 
consistent way. The SN-interface is assumed to be 
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ideal, and the outer boundary of the normal metal 
is specularly reflecting. The normal metal is char- 
acterized by a vanishing pairing interaction. We 
have neglected self-consistency of the pair poten- 
tial A(x) in the superconductor. The slight sup- 
pression of A(x) due to the proximity effect does 
not affect the low-energy spectrum of the normal 
metal since we assume d 3> Co = ^ v f /A. 

As a result of a detailed numerical study, we 
find that a gap forms at arbitrarily small impurity 
concentrations. This is shown in the upper panel 
of Fig. 1: even for values of the elastic mean free 
path I that are 50 times larger than the normal- 
layer thickness, the formation of the low-energy 
gap is clearly visible. The gap increases with l/l, 
saturates for I ~ d and then decreases again as 
expected from the dirty-limit theory since D ~ I. 
The gap does not depend on the location in the 
normal metal as can be seen in the lower panel 
of Fig. 1, i.e., it is a global feature. The shape of 
the LDOS, i.e., its dependence on energy, however, 
varies on traversing the normal layer. 

Figure 2 shows the dependence of the minigap 
on d/l. The behavior can be understood in a quali- 
tative way as follows: the linear dependence of the 
LDOS on energy in the clean case is caused by the 
existence of Andreev bound states with energies 
E n {6) = irhv F (n + l/2)cosfl/2d with 9 ~ tt/2. 
Here, 9 denotes the angle of a semiclassical trajec- 
tory with the normal vector of the interface. These 
trajectories (V-shaped structures bounded by two 
Andreev reflections and one specular reflection at 
the outer interface) have a length ~ d/ cos#. In al- 
most clean systems, this length is cut off by the 
mean free path, hence there is a lower energy scale 
of Ecap ~ Hvf/1 ~ Ta&/1 where Ta = Kvp/Zird 
is the Andreev temperature. The appropriate en- 
ergy scale for the diffusive limit can be obtained as 
h/t where t is defined by d <~ \[Dt. As a result, we 
obtain E Gap ~ fiD/d 2 ~ T A l/d. 

In conclusion, we have solved the real-time Eilcn- 
berger equation numerically and have determined 
the local density of states of a proximity sandwich. 
We find a hard gap that opens for an arbitrarily 
small concentration of impurities and is maximal 
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Fig. 1. Local density of states of the normal side of a 
proximity sandwich, d is the normal-layer thickness and 
I the elastic mean free path (both in N and in S); 
Ta = hvp /2nd) ; d = 1000£o- Upper panel: Adding impuri- 
ties leads to the formation of a minigap. Lower panel: The 
minigap is constant throughout the normal metal, but the 
energy dependence of the LDOS changes with location. 
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Fig. 2. Size of the minigap. Diamonds: numerical results. 
Dotted line: estimate in the almost clean case as given in 
the text. Dashed line: numerical fit to the dirty-limit results 
obtained in [4]. Inset: geometry of the SN-junction. 

if the elastic mean free path is of the order of the 
normal-layer thickness. 
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